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ABSTRACT 



We have used the IRAM Plateau de Bure Interferometer to map CO(3-2) 
emission from the gravitationally lensed Lyman break galaxy MS 1512-cB58. 
This is the first detection of a molecular emission line in any Lyman break system; 
its integrated intensity implies a total molecular gas mass of 6.6^| xlO 9 ho 2 7 M Q , 
while its width implies a dynamical mass of l-CTt^ x 10 10 csc 2 i \ M Q (for a 
flat f^A = 0.7 cosmology). These estimates are in excellent concordance with 
nearly all parameters of the system measured at other wavelengths, and yield a 
consistent picture of past and future star formation with no obvious discrepancies 
requiring explanation by differential lensing. In particular, we find that the age 
and remaining lifetime of the current episode of star formation are likely to 
be similar; the surface densities of star formation and molecular gas mass are 
related by a Schmidt law; and the fraction of baryonic mass already converted 
into stars is sufficient to account for the observed enrichment of the interstellar 
medium to 0.4 Z & . Barring substantial gas inflow or a major merger, the stars 
forming in the current episode will have mass and coevality at z = similar to 
those of a spiral bulge. Assuming cB58 is a typical Lyman break galaxy apart 
from its magnification, its global parameters suggest that the prescriptions for 
star formation used in some semi-analytic models of galaxy evolution require 
moderate revision, although the general prediction that gas mass fraction should 
increase with redshift is validated. The length of cB58's star formation episode 



1 Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by the INSU/CNRS (France), the MPG (Germany), and the IGN (Spain). 

2 Also: Department of Physics, 366 Le Conte Hall, University of California, Berkeley, CA 94720-7300 



- 2 - 



relative to the time elapsed over the redshift range 2.5 < z < 3.5 strongly argues 
against scenarios in which observed LBGs cohabit their halos with a large number 
of similar but "dormant" systems whose starbursts have faded or not yet begun. 
As a useful empirical result, we find that the observed line/continuum ratio 
for cB58 is similar to those of high-redshift systems with quite different dust 
luminosities and nuclear activity levels. Finally, we report the detection of a 
second source close to the position of the cD elliptical in the z = 0.37 lensing 
cluster, which may be nonthermal continuum emission from the cD or CO line 
emission from a hitherto unknown background galaxy at z ~ 1.48 or z ~ 2.73. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: kinematics 
and dynamics — galaxies: starburst — cosmology: observations 

1. Introduction 

The Lyman break technique (Steidel & Hamilton 1992) has proved to be extremely 
efficient in identifying large numbers of actively star-forming galaxies at z > 2 (Steidel et al. 
1999, 2003). While its selection criteria require only that Lyman break galaxies (LBGs) have 
rest-UV colors comparable to those of UV-bright starbursts at z = (e.g., Meurer, Heckman, 
& Calzetti 1999), the method has proved to identify high-redshift systems that resemble these 
local counterparts in other respects as well. The two populations have similar bolometric 
surface brightnesses (Meurer et al. 1997), UV-through-optical spectral energy distributions 
(Ellingson et al. 1996; Sawicki & Yee 1998), and blueshifted interstellar absorption from 
outflowing gas (Pettini et al. 1998, 2001). The primary difference is that LBGs appear to be 
forming stars at rates which are larger in proportion to their larger physical sizes (Meurer et 
al. 1997). Even the most careful studies of LBGs at relatively short wavelengths, however, 
provide only a partial picture of their important physical parameters. Although population 
synthesis models (Sawicki & Yee 1998; Papovich, Dickinson, & Ferguson 2001; Shapley et 
al. 2001) can reconstruct their past star formation histories, and rest-UV fluxes can be used 
to derive their present star formation rates (Madau et al. 1996; Steidel et al. 1999) modulo 
uncertain corrections for extinction (Ouchi et al. 1999; Peacock et al. 2000; Adelberger & 
Steidel 2000; Chapman et al. 2000; Seibert, Heckman, & Meurer 2002; Nandra et al. 2002; 
Webb et al. 2003), ignorance of their molecular gas contents makes it impossible to predict 
LBGs' future courses of star formation. Moreover, estimates of LBGs' total dynamical masses 
remain uncertain insofar as they are based solely on observations of nebular emission lines: 
ionized gas is more likely than molecular gas to reflect non-gravitational influences (Heckman 
et al. 1998), and less likely to trace the full rotation curve in local starbursts (Lehnert & 
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Heckman 1996). The detection of molecular gas in one or more LBGs is therefore of great 
interest in complementing the extensive work in the rest UV and optical. 

Within the LBG population, the brightest observable object is the gravitationally lensed 
system MS 1512-cB58 (hereafter cB58: Yee et al. 1996). Thanks to magnification of a 
background z = 2.7 source into a fold arc by the z = 0.37 cluster MS 1512+36, this galaxy 
has permitted detailed studies in the rest-frame UV (Ellingson et al. 1996; Pettini et al. 
2000, 2002) and optical (Bechtold et al. 1997; Teplitz et al. 2000), as well as detections in 
the rest-frame near (Bechtold et al. 1998) and far infrared (Baker et al. 2001; van der Werf 
et al. 2001; cf. Sawicki 2000). Two initial searches for CO (3-2) emission from cB58 were 
unsuccessful (Frayer et al. 1997; Nakanishi et al. 1997); notwithstanding cB58's sub-solar 
metallicity, these raised the question of how stars could be forming at such a seemingly 
high rate out of such a seemingly small reservoir of molecular gas. Our measurement of a 
surprisingly low 322 /im flux density (A Q bs = 1200 /im: Baker et al. 2001), implying a lower 
UV extinction and total star formation rate than those assumed by Frayer et al. (1997), 
suggested that molecular line emission might also be weaker than previously expected. The 
measurement by Teplitz et al. (2000) of nebular emission at a redshift apparently outside 
the bandpass of the Nakanishi et al. (1997) observations strengthened our conviction that a 
new and more sensitive search for molecular gas in cB58 was warranted. 

In this paper, we report our successful detection of CO (3-2) emission from cB58- in 
fact, the first detection of molecular line emission from any LBG. §2 describes our acquisition 
and reduction of the data. §3 discusses our measurement of the CO (3-2) line parameters for 
cB58 and our derivations of molecular gas and dynamical masses. §4 presents our detection 
of an additional source "X" in the field, close to the position of the cD elliptical in the 
z = 0.37 lensing cluster. We consider the likeliest origin of this emission to be nonthermal 
continuum emission from the cD or CO line emission from a background galaxy, although 
neither explanation is entirely satisfactory. §5 lays out a coherent picture of the past, present, 
and future of star formation in cB58, discusses the implications for our understanding of 
the z ~ 3 Lyman break population as a whole, and assesses source-to-source variation 
in the millimeter line/continuum ratio at high redshift. §6 summarizes our conclusions. 
Throughout this paper, we will assume H = 70 /10.7 km s _1 Mpc -1 and a flat Q\ = 0.7 
cosmology. 3 The correction of observed fluxes for lensing magnification generally makes use 
of the most satisfactory model MVIa of Seitz et al. (1998 = S98). In this interpretation, all 
of the background rest-UV source is magnified by a factor 2.39 into a counterimage (A2), 
while a lesser fraction (/ arc ~ 0.54 — 0.67) of the background source is magnified by a larger 



3 For cB58 at z — 2.729, this cosmology gives £>l = 22.7/i *Gpc and Da — 1.63/i ?Gpc. For the lensing 
cluster at z = 0.37, D L = 1.98 Gpc and D A = 1.06 /iq? Gpc. 
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factor (oc f aT l ~ 59 — 48) into the cB58 fold arc, giving an arc/intrinsic brightness ratio of 
M VY - 31.8. 

2. Observations 

We observed cB58 with the IRAM Plateau de Bure Interferometer (PdBI: Guilloteau et 
al. 1992) for 18 intervals from April 2001 through August 2002. The array included between 
four and six 15 m diameter antennas; during our observations, these were arranged in the 
compact D configuration whose 15 baselines range from 24 to 82 m in length. Each antenna 
is equipped with both 3 mm and 1 mm SIS receivers that can be used simultaneously; these 
were tuned in single-sideband mode at 3 mm for 11 of the 18 observing epochs (encom- 
passing ~ 60% of the data), and in double-sideband mode for the remaining 3mm and 
all 1mm data. Receiver temperatures were 60-80 K, and corrected outside the atmosphere 
yielded typical system temperatures of ~ 120 K and ~ 400 K at 3 mm and 1 mm, respec- 
tively. For our observations of cB58, we adopted a phase center (aj2ooo = 15 h 14 m 22.25 s and 
5 J2000 = +36°36'24.40") 0.6" east of the coordinates listed in NED (Abraham et al. 1998). 
For the z H n = 2.7290 ± 0.0007 reported by Teplitz et al. (2000), the CO(3-2) and CO(7-6) 
rotational transitions are redshifted to 92.73 GHz and 216.32 GHz. We deployed four cor- 
relator modules- giving a total of 580 MHz of continuous bandwidth at 2.5 MHz channel -1 
resolution- across each of the two lines; both were observed in the lower sideband relative 
to the reference frequency. After reducing the 2001 subset of the data, we shifted the 3 mm 
receiver tuning to 92.83 GHz, in order to obtain more symmetric coverage of the already 
apparent CO(3-2) line at z C q < -^hii (see §3.1). As a result, only the central ~ 480 MHz of 
3 mm bandwidth containing data from all 18 observing epochs has been used to produce the 
line and continuum maps discussed below. 

We calibrated the data using the CLIC routines in the IRAM GILDAS package (Guil- 
loteau & Lucas 2000). Phase and amplitude variations within each track were calibrated out 
by interleaving observations of 3C345 (at a separation A9 = 17.7° from cB58 on the sky) 
every 30 minutes; we also observed the weaker but closer (A9 = 1.9°) quasar 1504+377 to 
permit a check of our astrometric uncertainty (see §3.1 below). Passband calibration used 
3C345 or another bright quasar. The overall flux scale for each epoch was set by comparing 
3C345 with 3C273 and/or with model sources CRL618 and MWC349. Reference fluxes for 
3C273 were extracted from an archive of dedicated monitoring data from the PdBI and the 
IRAM 30 m telescope. We estimate the accuracy of our final flux scales to be ~ 15%. 

Before constructing uv tables from our calibrated data, we smoothed them to a frequency 
resolution of 5 MHz (i.e., velocity resolutions 16.2 km s" 1 for the CO(3-2) line and 6.9kms _1 
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for the CO (7-6) line). After editing the data for quality in the DIFMAP package (Shepherd 
1997), we were left with the (on-source, six-telescope array) equivalents of 32 hours of 
3 mm data and 25 hours of 1 mm data. We inverted the data cubes with natural weighting 
(but no deconvolution) using the IMAGR task in the NRAO AIPS package (van Moorsel, 
Kemball, & Greisen 1996). At 3 mm, the lower-sideband data yield a synthesized beam of 
8.2" x 4.8" at P. A. 73°, and a noise level per 16.2kms _1 channel that increases with velocity 
(decreases with frequency) from 0.8 to 1.1 mJy beam -1 across the 480 MHz bandpass. The 
upper-sideband 3 mm data from the epochs with double-sideband tuning give a synthesized 
beam of 8.2" x 4.6" at P. A. 75°, and a sensitivity to 95.7 GHz continuum emission of 
0.24 mJy beam -1 . At 1mm, the lower-sideband data give a synthesized beam of 3.0" x 2.1" 
at P. A. 91°, and a noise level per 6.9 km s -1 channel of 5.3 mJy beam -1 across the entire 
580 MHz bandpass. When both 1mm sidebands are averaged together, the sensitivity to 
217.8 GHz continuum emission is 0.41 mJybeam -1 for the same synthesized beam. 



3. CO(3-2) emission from cB58 

3.1. Line parameters 

Cursory inspection of the 3 mm data cube reveals the presence of a line emission source 
close to the phase center with velocity width ~ 200 km s" 1 . To quantify the significance and 
parameters of this source, we convolved the 16.2 km s" 1 resolution cube with a twelve-channel 
boxcar filter. The resulting cube contains 194.0 km s™ 1 channels, each of which overlaps its 
immediate neighbors by 177.8 km s _1 . The strongest source anywhere in this smoothed cube 
is shown in Figure 1: neglecting primary beam correction, it is a 1.7 ± 0.3 mJy beam -1 
peak centered at 92.79GHz (i.e., ~ 190 kms -1 blueward of the nominal znu). A two- 
dimensional Gaussian fit in the image plane indicates that it is spatially unresolved; Table 1 
lists the position derived from a fit that fixes the dimensions of the Gaussian to those of the 
synthesized beam. Two main lines of evidence argue that this 1.7 ± 0.3 mJy beam -1 peak is 
in fact the 0.33 ± 0.06 Jykms -1 core of a CO (3-2) emission line from cB58: 

1. The position of the source. We can accurately register optical source positions into the 
radio reference frame for this field by identifying the (17er) 7mJy radio core in the 0.4" 
resolution 6.1cm map of Laurent-Muehleisen et al. (1997) with the central cD galaxy 
in the optical catalogue of Abraham et al. (1998). 4 Under this assumption, the position 



4 NED positions for sources in this field have been calculated by identifying the cD with a less precisely 
located 21cm source (Becker, White, & Helfand 1995). 
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of the optical fold arc agrees with the position of the putative CO (3-2) detection to 
within the latter's astrometric uncertainties (see Table 1). Following Downes et al. 
(1999), we have calculated these uncertainties as the sum of three terms in quadrature: 

(a) The ~ 0.01" absolute uncertainty of the PdBI reference frame; this stems from 
the fundamental uncertainties of quasar positions. 

(b) The relative uncertainty introduced by using a calibrator (3C345) that is offset by 
A6 ~ 17.7° from cB58. In principle, this uncertainty is about 0.01" per degree of 
separation, i.e., ~ 0.17". However, because any error is expected to vary smoothly 
with position, we can exploit the fact that the quasar 1504+377 (only 1.9° away 
from cB58, but 19.0° away from 3C345) lies within 0.125" of the phase center in 
a map also calibrated by 3C345. For cB58, lying between 3C345 and 1504+377, 
we can then impose the somewhat tighter limit of (17.7/19.0) x 0.125" ~ 0.12". 

(c) The accuracy of our centroiding, which AIPS indicates is ±0.6" in right ascension 
and ±0.4" in declination. This term dominates the other two. 

2. The redshift of the line. Pettini et al. (2000, 2002) have arrived at estimates of the 
systemic redshift of cB58 from rest-UV photospheric absorption features that are blue- 
ward of the Teplitz et al. (2000) Zhh by ~ 110 — 180 km s -1 in the rest frame of the 
source (see Table 2). From fits to the CO line profile (see below), we estimate a 
zco = 2.7265^q;qqq5> whose ~ 200 km s" 1 blueshift relative to ^hii agrees fairly well 
with the rest-UV results. We do not view the consistent difference from zun as nec- 
essarily astrophysically significant; if Ha were to arise largely in a wind, for example, 
we would expect zun to be lower than zqo, contrary to what is seen. The discrepancy 
may rather reflect simple uncertainties in wavelength calibration. 



Three further consistency checks offer no evidence against this proposed identification. First, 
our detection is consistent with the Frayer et al. (1997) 3<r upper limit of 0.86 Jykms -1 for 
a 194.0 km s _1 velocity width. Second, we see no sign of a velocity gradient in the CO(3-2) 
data cube; this negative result is consistent with the absence of any shear along the arc 
in rest-UV observations of higher spatial resolution (Pettini et al. 2002). Finally, the lack 
of obvious 1mm continuum emission at the position of the putative CO (3-2) line (above 
a 3(7 upper limit of 1.2 mJy) agrees with the results of Baker et al. (2001), who report a 
250 GHz continuum flux density for cB58 of 1.06 mJy. Adopting their suggested T d = 33 K 
dust temperature and e u = 2 dust emissivity index, the 217.8 GHz flux density expected in 
our PdBI map would be only 0.68 mJy (see also §4.1). 

Given this identification, we proceed to calculate the parameters of the CO (3-2) emis- 
sion in cB58. While Figure 1 gives a fair estimate of the flux in the core of the line 



7 



(~ 0.33 Jykms -1 ), it fails to account for additional emission in the line wings extending 
blueward and redward of this particular 194.0 km s -1 chunk of the spectrum. We have there- 
fore fit Gaussians to the line profile through the brightest pixel in Figure 1, as extracted from 
versions of the original 16.2 km s -1 cube that have been rebinned to 48.5 km s -1 resolution 
and corrected for primary beam attenuation (e.g., Figure 2). To estimate the uncertainties in 
the line widths and intensities, we have performed these fits on spectra smoothed from 16.2 
to 48.5 km s -1 with all three possible choices of rebinning (i.e., averaging channels {2-4, 5-7, 
...} and {3-5, 6-8, ...} as well as {1-3, 4-6, ...}). In addition to z C q = 2.7265, we estimate a 
total line flux F CO (3-2) = 0.37±0. 08 Jykms -1 and a velocity dispersion aco = 74±18kms -1 
(equivalently, FWHM co = 174 ± 43 km s -1 ). 

In the 1mm data cube, we have summed putative CO(7-6) channels over a 194.0 km s -1 
velocity interval within 2.3 km s -1 of that integrated to produce Figure 1. At the position 
of the CO (3-2) peak, no CO (7-6) emission is seen; the nominal line flux is Fco(7-6) = 
0.15 ±0.21 Jykms -1 . For a fair estimate of the CO (7-6) /CO (3-2) integrated intensity ratio 
7^7/3, however, we must allow for the possibility that relative weakness in the CO (7-6) line 
can partly result from a relative lack of uv sampling at short r uv = B/\ b s , B denoting 
the projected baseline length. We have therefore truncated both of the original CO (3-2) 
and CO(7-6) uv datasets below an identical minimum r uv — 1.1 x 10 4 , and reinverted them 
with natural weighting. After correcting for primary beam attenuation and convolving the 
new CO(7-6) map to the 7.2" x 4.1" (P. A. 80.9°) beam of the new CO(3-2) map, we take 
their quotient. In flux units, the line ratio in a small region centered on the position of the 
CO(3-2) line detection is then 0.16 ± 0.04. Translated to brightness temperature units, this 
yields a nominal integrated intensity ratio 



The lack of CO(7-6) emission in a noisy 1mm map means that this value may not be 
a particularly meaningful upper limit. However, we can safely conclude that there is no 
evidence for large amounts of warm, dense molecular gas in cB58 like those observed in the 
submillimeter galaxy SMM J14011+0252 at a similar redshift (Downes & Solomon 2003). 

To evaluate whether the width and intensity of the CO (3-2) line in cB58 are consistent 
with what we would expect for a relatively unobscured starburst, we must compare it to the 
parameters of CO (3-2) emission observed in appropriate z = analogues. Mauersberger 
et al. (1999) report single-pointing CO(3-2) observations of a wide range of galaxy types; 
however, a better comparison sample is that of Dumke et al. (2001), who have mapped the 
CO (3-2) line across the disks of twelve nearby starburst and star-forming galaxies. The most 
revealing quantity for comparison with cB58 is the beam-diluted brightness temperature 




(1) 
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/A?b, where fx is the areal filling factor of a galaxy's molecular clouds and T b is the rest- 
frame brightness temperature of a single cloud. Following Downes & Solomon (2003), we can 
write the "unlensed" CO luminosity in terms of both observed and rest-frame quantities: 

L' co _/ Fen \ /A ohs \ 2 / D A \2 ( Mt 
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(jyl^J (m^J (mp^J (—J + (2) 
7b \ / At; \ /fy\ /fA\ t Attco \ ( Da \ 2 (3) 
KJ Vkms- 1 / U/lli Varcsec 2 / vMpc/ 1 ' 



Kkms _1 pc 2 Vjykms" 1 / Vmm/ VMpc 

= 23.5 



where Af2 co is the solid angle in the source plane, / v is a velocity filling factor generally 
assumed to be unity due to the smoothness of extragalactic CO line profiles (e.g., Dickman, 
Snell, & Schloerb 1986), and .Mco is the factor by which the intrinsic CO(3-2) emission is 
magnified due to lensing. Solving for f\T h , we find 

f A T h ( F co \f\ oh s\ 2 f Av x-y A^co A-V-McoA- 1 ,, , , (A , 

IT = 15 -H J^knT^) vn^J IknT^J la^J (— ) (1 + z) (4) 

independent of cosmology. For the Dumke et al. (2001) sample, we can set Mco = 1 + 2 = 1 
and A bs = 0.87mm; the CO(3-2) line flux, line width, and solid angle 5 of the dominant 
"central peak" they determine for each system are listed in their Table 3. The resulting 
{/a 7b} range from 0.4 to 5.4 K. For cB58, we can assume that the regions producing CO(3- 
2) line and rest-UV continuum emission (i.e., molecular clouds and the stars forming from 
them) are very similar, so that M.co — M-uv — 31.8. By the same token, since model 
MVIa of S98 implies that the background LBG lensed into the cB58 arc has semi-major 
axes a uv x &uv — 0.25" x (0.25/2.39)", we can assume AQ co ~ Af2 uv = 0.082 arcsec 2 . For 
A bs = 3.23 mm and the line flux and FWHM derived above, Equation 4 then implies cB58 
has f A T h ~ 0.5lg;|K. This value falls near those seen in the centers of M51 (0.4), NGC278 
(0.6), and NGC4631 (0.8) in the Dumke et al. (2001) sample. 

The agreement of f\Th estimates for cB58 and three nearby star-forming galaxies is 
strictly empirical- in particular, we have not (yet) invoked any relationship between L' co ^_ 2 ^ 
and molecular gas mass- and has three general implications. First, the consistency gives 
tentative support to our assumption that the system's rest-UV continuum and CO line 
emission are identically lensed; there is no order-of-magnitude discrepancy that can only be 
accounted for by differential lensing. Second, it appears entirely plausible that cB58 contains 
molecular clouds whose distribution (i.e., /a) an d typical internal conditions (i.e., Tb) overlap 
with the ranges seen in local star-forming galaxy populations. A stronger formulation of 
this point is not warranted, given that the relative weighting of f\ and Tb in the product 



5 Computed as ^D\Dij\ from source diameters D\ and D2 in arcseconds. 
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of both could still differ at high and low redshift. Finally, the molecular ISM in cB58 does 
not appear to resemble that of an ultra-luminous infrared galaxy (ULIRG) merger. Five of 
the six ULIRG molecular disks from which Lisenfeld et al. (1996) detect CO (3-2) emission 
have CO(l-O) line widths measured by Solomon et al. (1997) and intrinsic CO(l-O) and/or 
CO (2-1) source sizes inferred by Downes & Solomon (1998). Assuming that these systems' 
Av and Af2co are similar for all CO transitions, we estimate values for their f\Tb in the 
range 16 — 47 K. By the same token, the parameters of CO (3-2) emission in the lensed IR- 
luminous galaxy F10214+4724 (z = 2.3) mapped by Downes, Solomon, & Radford (1995) 
imply fxTb — 10 K- a factor of 20 larger than in cB58. 



The CO(l-O) line luminosity is well established as a proxy for the total gas mass grav- 
itationally bound in a galaxy's molecular clouds (Dickman et al. 1986; Solomon & Bar- 
rett 1991; Wolfire, Hollenbach, & Tielens 1993). In terms of a conversion factor X co = 
NnJ I C o(i-o) cm" 2 (Kkms" 1 )" 1 , 



after including a factor of 1.36 to account for helium. Although Xco — (2 — 3) x 10 20 
applies reasonably well for molecular clouds in the Galaxy and in the disks of nearby spirals 
(Strong et al. 1988; Solomon & Barrett 1991; Scoville & Young 1991), for galaxies with 
subsolar metallicity the universality of the conversion factor is a contested issue. Wilson 
(1995) and Arimoto, Sofue, & Tsujimoto (1996) argue that Xco rises for lower-metallicity 
systems; averaging their recommended values for the metallicity of cB58, which we take 
to be 12 + log(0/H) ~ 8.39 from Teplitz et al. (2000), would imply X co = 7.3 x 10 20 . 
In contrast, Rosolowsky et al. (2003) find no systematic trend of X C o with metallicity in 
the range 0.25 < Z/Z & < 1 across the disk of M33, while Bolatto et al. (2003) report 
X co = 7.0 x 10 20 for the Z/Z G = 0.1 SMC, considerably smaller than the Wilson (1995) 
and Arimoto et al. (1996) relations would predict. To reflect the range of these results, we 
adopt X C o = (4.5 ± 2.5) x 10 20 cm~ 2 (Kkms -1 )" 1 for cB58 in this paper. 

Next, we must estimate cB58's integrated intensity ratio 72-3/1 = 7co(3-2)/7co(i-o)- 
Dumke et al. (2001) report values ranging from 0.4 to 1.4 in the centers, and from 0.5 
to 1.0 in the disks, of their sample of twelve galaxies. The centers of M51 and NGC278 are 
observed to have 72 3 /i = 0.5 and 0.8, respectively; based on the analogy with cB58 in §3.1 



3.2. Molecular gas mass 




above, we adopt IZ3/1 = 0.65 for cB58. Now combining Equations 2 and 5, 




(6) 



Again assuming A^co — M-vv-, and including the uncertainty in the appropriate value of 
Xqo, we estimate for cB58 M gas = 6. 6*4 3 x 10 9 M@. We make no correction for the 
atomic mass component of cB58's ISM, although we note that in low-metallicity systems, a 
large fraction of the gas that is gravitationally bound into the "molecular" clouds traced by 
the CO lines (and thus readily available for star formation) may already be atomic in phase. 

cB58's molecular gas mass is comparable to those seen in the more gas-rich members 
of the local spiral population (e.g., Heifer et al. 2003), somewhat smaller than those seen in 
nearby ULIRGs (Bryant & Scoville 1999; cf. Downes & Solomon 1998), and up to an order 
of magnitude less than those inferred for submillimeter-selected systems at comparably high 
z > 2.5 (Frayer et al. 1998, 1999; Neri et al. 2003). If its molecular gas is spread over a disk of 
radius r co = 0.25" ~ 2.0 /i^kpc in the source plane (see §3.3 below), the corresponding mass 
surface density will be S gas = 540^0 M© pc -2 . Such surface densities are more reminiscent 
of the centers of gas-rich spirals (Heifer et al. 2003) than of ULIRGs, consistent with the 
similarity noted above between cB58 and the Dumke et al. (2001) sample in terms of f^Th- 
Since the use of a single conversion factor Xqq in different galaxies relies on the assumption 
that their molecular cloud properties are reasonably similar (Dickman et al. 1986), this 
consistency serves to validate our use of a (metallicity-corrected) Galactic Xco to estimate 
M gas in the first place. 6 



Our spectrally resolved CO (3-2) detection means that we can estimate the dynamical 
mass of cB58 from its velocity width according to 



Here a is a constant of proportionality whose value depends on the intrinsic geometry of 
cB58's molecular gas, and on how that geometry couples with the effects of lensing. Since 



6 Indecd, our use of a conversion factor exceeding the Galactic Xqo by a factor of 2.25 implies that 
in terms of its mean molecular hydrogen column (^Vh 2 °^ -^Co/A^tOi cB58 will more closely resemble the 
galaxies in the Dumke et al. (2001) sample whose beam-diluted brightness temperatures are more typical 
(i.e., not among the lowest observed). 



3.3. Dynamical mass 
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molecular gas is extremely dissipative and unlikely to remain molecular in phase if severely 
disturbed, it is more likely than not that the CO emission from cB58 arises in a rotating disk. 
Indeed, the S98 source plane reconstruction suggests that the fraction (/ arc ~ 0.6) of the 
background LBG being magnified into the fold arc corresponds to half of such a disk, slightly 
inclined and with half-light radius ruv — 0.25". In this scenario, assuming that rco — r uv 
and that the magnified half of the disk traces only half of the galaxy's full line width, we would 
infer a ~ 4csc 2 i. This corresponds to a dynamical mass M dyn = 1.01^4 x W 10 csc 2 i h^} M Q 
within a radius of 2.0 kpc, with the uncertainties in M dyn reflecting only the uncertainties 
in oco 

If we wish to avoid relying too heavily on the details of the source morphology re- 
constructed by S98, we can also consider the possibility that cB58's molecular gas has a 
spherical rather than disklike distribution. In this scenario, a = 5 (for a sphere of uniform 
density) and correction of the observed half-light radius of A2 for lensing yields a character- 
istic r co — (0.25/V2.39)" = 0.16". The implied dynamical mass is now a somewhat smaller 
M Ayn = 0.8±° i x lO lo /i a ^M within a radius of 1.3 kpc. We view this estimate as less 
plausible than the result for a disklike geometry, however, for the a priori reasons noted 
above. 

For either assumed gas configuration, we can compare Md yn with the M gas derived in 
§3.2 above to estimate the molecular gas mass fraction / gas in the center of cB58. For a 
disk geometry, / gas = 0.66^q'| 2 sin 2 i \ within a 2.0 h$ \ kpc radius; for a spherical geometry, 
/ gas = 0.82+°:64^a7 within 1.3 kpc. Both scenarios imply a rather high gas mass fraction, 
although the uncertainties in both estimates (including sin 2 i for the disk case) make it 
impossible to exclude the possibility that / gas <C 0.5. 

4. Emission from a second source 

4.1. X and its possible origins 

Further examination of our 3 mm data cube- specifically, the velocity channels outside 
the range spanned by cB58's CO (3-2) emission- reveals a second source in the same field. 
Averaging over a 200 MHz bandwidth centered at 92.96 GHz (i.e., over a 645 km s -1 velocity 
range lying 523 km s -1 blueward of zco for cB58) yields the map shown in the center panel 
of Figure 3. The new source (hereafter "X") has an apparent flux density 0.87 ± 0.16 mJy 
(±0.13 mJy for the uncertain flux scale) after primary beam correction. A Gaussian fit 
gives the position listed in Table 1, which coincides with the radio continuum core of the 
cD in the lensing cluster (Becker et al. 1995; Laurent-Muehleisen et al. 1997) and is clearly 
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distinct from the position of cB58's CO(3-2) centroid. 7 Puzzlingly, however, we do not detect 
comparably significant emission from X in a 140 MHz bandwidth centered at 92.66 GHz 
(nominally 0.46 ± 0.24 (±0.07) mJy), or in our 95.7 GHz continuum map (nominally 0.26 ± 
0.24 (±0.04) mJy), as shown in the left and right panels of Figure 3. X also fails to appear 
in our 217.8 GHz continuum map (nominally 0.79 ± 0.45 (±0.12) mJy). 

Since the noise is 50% higher in our maps at 92.66 GHz (due to a smaller bandwidth and 
the noise gradient described in §2) and at 95.7 GHz (due to the dearth of upper-sideband 
data) than at 92.96 GHz, it remains possible that X is indeed a weak continuum source. In 
this case, its positional coincidence with the cD suggests a likely nonthermal origin. The 
cD's radio core has a 1.4 GHz flux density ~ 1 mJy at 5" resolution (Becker et al. 1995), 
and a 4.9 GHz flux density ~ 7mJy at 0.4" resolution (Laurent-Muehleisen et al. 1997), 
meaning it either varies or has an unusually inverted slope (cxf^ ~ —1.56 for S u oc v~ a ). In 
the Dallacasa et al. (2000) sample of High Frequency Peakers (HFPs), reaches values 
as extreme as —1.72, so the latter possibility is plausible. Assuming the cD radio core 
is young, with a ~ —1.56 below a turnover frequency z/ pea k as high as those of the most 
extreme HFPs (~ 22 GHz) and a ~ ±0.8 above z/ pea k (e.g., Condon 1992), we would predict 
a 92.96 GHz flux density ~ 23mJy. A less aggressively tuned SED could still account for 
all of the observed 0.87 mJy from X. This scenario would explain why no emission from X 
appears in our 217.8 GHz continuum map, since a spectral index a = ±0.8 would imply 
5*217.8 — 0.44 mJy, i.e., at the level of the noise. 

A less plausible continuum scenario accounts for X as thermal dust emission associated 
with star formation in the cD. Analogous systems do exist; of the two 850 /im detections of 
z ~ 0.25 cDs by Edge et al. (1999), for example, at least one (in Abell 1835) has an unam- 
bigous dust origin. However, the tight constraints on emission from X at i/ > 92.96 GHz 
make it unlikely that the cD in MS 1512±36 is similar. For the median dust temperature 
T d = 58 K and (relatively flat) emissivity index e u = 1.35 in the local starburst sample of 
Yun & Carilli (2002), a z = 0.37 galaxy with flux density 0.87mJy at 92.96 GHz should 
have flux densities ~ 14, 22, and 65mJy at 1.4mm (217.8GHz), 1.2mm, and 850 pm, re- 
spectively. It is marginally possible that our PdBI map could have resolved out 14mJy of 
217.8 GHz emission from X. It is completely inconceivable, however, that bolometer array 
observations of cB58 could have missed such strong emission from a second source so close 
on the sky. For the 11" and 15" (HPBW) beams of MAMBO at the IRAM 30 m telescope 
and SCUBA at the JCMT, a 4.8" separation translates to attenuation by factors of 0.59 and 



7 The tight astrometric constraints in a common radio reference frame also imply that cB58 cannot be 
the source of both the CO(3-2) emission discussed in §3 and the 6.1cm continuum emission mapped by 
Laurcnt-Muchlciscn et al. (1997). 
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0.75, respectively. Accounting for this correction, Baker et al. (2001) should have measured 
~ 13 mJy at 1.2 mm for cB58, rather than the 1.06 ± 0.35 mJy observed, while van der Werf 
et al. (2001) and Sawicki (2001) should have measured ~ 49 mJy at 850 /xm, rather than their 
4.2 ± 0.9 mJy detection and 3.9 mJy 3<x upper limit. Similar inconsistencies are encountered 
for any reasonable choices of and e v . 

If we instead take Figure 3 at face value and interpret X as an emission-line source, 
its positional coincidence with the cD must be essentially accidental. Its irregular velocity 
"profile" (Figure 4) does not permit a detailed fit for an exact centroid or width, although it 
is clear that the latter must be comparable to the bandwidth smoothed to obtain the center 
panel of Figure 3. Such a ~ 645 km s" 1 velocity width is incompatible with the 246kms _1 
velocity dispersion for the cD required by the S98 lensing model; moreover, a rest frequency 
~ 127.35 GHz does not correspond to any strong recombination line (Towle, Feldman, & 
Watson 1996) or molecular rotational transition. We therefore adopt the reverse strategy 
of assuming X is one of the strongest extragalactic emission lines- i.e., a CO rotational 
transition- and asking whether an object at the corresponding redshift can provide a self- 
consistent explanation. As CO(1-0) emission, X would necessarily come from a z ~ 0.24 
galaxy lying in front of the cluster. The only possible candidate is the foreground spiral "S" 
seen ~ 4.4" southeast of X, which S98 suggest has 0.25 < zs < 0.37 on the basis of its angular 
size and V — R color. However, its apparently face-on orientation makes a ~ 645 km s" 1 
linewidth highly implausible, particularly at a radius ~ 16.6/iQ^kpc from the center of a 
spiral that would otherwise appear to be CO-poor. 

As CO(2-l) emission from a zx ~ 1-48 galaxy or CO(3-2) emission from a (further) 
zx ~ 2.73 galaxy, X would lie behind the lensing cluster and therefore undergo some degree 
of magnification. A significant challenge for this scenario is that any background source 
producing an image within an arcsecond of the central cD should also produce a number 
of counterimages (over the same velocity range) elsewhere in the cluster field. Velocity 
considerations alone mean X cannot be CO (3-2) emission from a fifth image in the W set 
of S98, since here the background source redshift lies only 150 km s -1 from that of cB58 
(Teplitz, Malkan & McLean 2002). To explore the possibility that X is CO(2-l) or CO(3-2) 
emission from a hitherto unknown background source, we have implemented the S98 lensing 
model within the gravlens package (Keeton 2001). As discussed in Appendix A, we have 
then assumed that X is one image of a lensed background point source and calculated where 
the corresponding counterimages should lie for both z x = 1-48 and zx = 2.73. The only 
configurations that do not predict a strong counterimage where none is observed require that 
X correspond to two merging images of opposite parity on opposite sides of a critical curve. 
Because the details of such models depend sensitively on the exact shape of the lens mass 
profile, we view our particular results as largely speculative. We do consider zx — 2.73 a more 
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appealing possibility than zx — 1.48, since we already know that cB58 and the background 
source producing the W image set lie at about the same redshift (Teplitz et al. 2002). The 
projected separation between the background galaxies lensed to produce cB58 and X would 
be only 3.1" ~ 25kpc in this case; the presence of an additional, more obscured and gas- 
rich companion that could have helped trigger cB58's star formation recalls the multiple- 
component morphology of the LBG Westphal-MMDll (Chapman et al. 2002). However, for 
either zx, a lensed background source would need to have sufficiently strong dust extinction 
to avoid easy detection behind the cD in HST optical imaging, yet sufficiently weak dust 
emission to avoid easy detection in 1.2 mm and 850 /xm bolometry. This combination of 
properties requires a significant, perhaps fatal degree of fine tuning. 



4.2. Implications of X for parameters of cB58 

There are two ways in which the presence of source X within 4.8" of cB58 on the 
sky might lead us to adjust the properties of the emission we attribute to the LBG. The 
first is the prospect that continuum (or broad line) emission from X has contaminated our 
measurements of the CO (3-2) emission line within the PdBI data cube. This turns out 
not to be a problem. If we treat X as line emission, then the best estimate of its velocity 
centroid and width (see Figure 4) place it mostly blueward of cB58's CO(3-2) emission, and 
the positional offset suppresses any remaining contamination. If we treat X as nonthermal 
continuum from the cD and subtract a model of the line-free channels from the uv data 
before inverting a now "line-only" cube, we measure CO (3-2) line parameters for cB58 that 
are virtually identical to those derived in §3 without the benefit of "continuum subtraction" 
(e.g., the peak in Figure 1 nominally increases from 1.7 ± 0.3 mJy to 1.8 ± 0.3 mJy). 

The second possible problem is that if X is continuum emission from the cD or line 
emission from a dusty background source, there may be some contamination of the 1.2 mm 
and 850 /im dust emission attributed to cB58 by Baker et al. (2001) and van der Werf et al. 
(2001), respectively. The amount of dust that could be associated with a lensed line source 
is difficult to quantify. The possible contamination by nonthermal continuum from the cD is 
more tractable: assuming a = +0.8 and accounting for the different MAMBO and SCUBA 
beam sizes, we would estimate that roughly 0.22 mJy of each detection could have come 
from the cD. This would represent ~ 5% of the 850 /im detection, but ~ 22% of the 1.2 mm 
flux density. To the extent that this has occurred, the weaker dust emission attributed to 
cB58 would support the conclusion of Sawicki (2001) that the galaxy's long-wavelength SED 
differs from those of local starbursts, although this difference would still be less dramatic 
than past (over)estimates of its rest-UV reddening have suggested (Baker et al. 2001). 
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5. Discussion 

5.1. Star formation in cB58 

Previous studies of cB58's rest-UV spectroscopic properties have placed strong con- 
straints on its most likely star formation history. The underabundance of nitrogen and 
iron-peak elements inferred from a Keck/ESI spectrum (Pettini et al. 2002; cf. Teplitz et al. 
2000) requires that star formation must have begun at some tsF(past) < 300 Myr, consistent 
with the lack of evidence for any older stellar population in rest-frame near-IR emission 
(Bechtold et al. 1998). The initial population synthesis analysis by Pettini et al. (2000) of 
their Keck/LRIS spectrum has been revisited by De Mello, Leitherer, & Heckman (2000) 
and Leitherer et al. (2001) using successive improvements of the original Starburst9 9 package 
(Leitherer et al. 1999). Overall, these three exercises find satisfactory agreement with the 
observed spectrum if cB58 has been forming stars at a constant rate with a Salpeter (1955) 
Initial Mass Function (IMF) up to 100 M for the last t S F( P ast) = 20 - 100 Myr. We can use 
this picture of cB58's recent past to estimate its star formation rate in the present. The 
Z = 0.4 Z Starburst99 models for the above range in £sF( P ast) predict 
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where we have extended the Salpeter IMF to 0.1 M Q . For cB58, Teplitz et al. (2000) measure 
a Balmer decrement of 3.23 and an Ha flux of (12.56 ± 0.37) x 10~ 16 erg cm -2 s -1 ; correcting 
the latter by a factor of 1.32 for LMC extinction (Howarth 1983) and assuming Ain a — 31.8 
implies SFR = (23.9 ±0.7) 2 M yr -1 . Given that the distributions of UV and Ha emission 
in local starbursts are broadly similar (Conselice et al. 2000), Aiua — M-vv should indeed 
hold if the fraction of the background source being lensed is truly as high as / arc ~ 0.6. 
The rate of star formation per unit area in cB58, Esfr = 1-9 ±0.1 M yr _1 kpc -2 , appears 
within the uncertainties to be following a Schmidt (1959) law dependence on S gas . From the 
best-fit relation determined by Kennicutt (1998), 

SsFR = 2.5 x 10- 4 (-^) 1A (9) 



MQvr-ikpc" 2 ' VM pc- 2 , 
we would have predicted cB58 to have Ssfr = 1.6± 2 ;2 M yr -1 kpc~ 2 . 

After continuous star formation at the current rate for the last 20-100 Myr, cB58 should 
already have produced M stars = (0.5 — 2.4) x 10 9 M Q in the present episode. For any such 
stellar mass, the implied baryonic mass- M bary = M gas ±M stars ~ (7.1 — 8.9) x 10 9 M - is still 
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safely less than the best estimate of the dynamical mass (M dyn = 1.0 x 10 10 csc 2 i Iiq} M ), an 
encouraging consistency check. For several reasons, we consider a stellar mass (oc isF(past)) & t 
the upper end of the allowed range to be the most plausible scenario. First, estimates of small 
tsF(past) (and low M stars ) that involve rest-UV continuum colors (e.g., Ellingson et al. 1996; 
Pettini et al. 2000) entail higher levels of extinction than are indicated by several other lines 
of evidence (Baker et al. 2001). Second, the difference Md yn — Mb ary should be small, since 
(a) there is no evidence for any very old stellar population (see above), and (b) we do not 
expect large dark matter contributions to the total mass budget at R < 2.0 h^kpc. Indeed, 
assuming % = 90° and tsF(past) — 140 Myr (all but indistinguishable from tsF(past) = 100 Myr 
in a population synthesis model with continuous star formation) would leave M dyn = M bary 
exactly. We will therefore assume tsF(past) — 140 Myr and M sta rs — 3.4 x 10 9 M in what 
follows. Third, if cB58 is a typical LBG in respects other than its fortuitous location behind a 
lensing cluster, we would expect its M stars to be comparable to those inferred from population 
synthesis analyses of LBGs with similar intrinsic brightnesses. Correcting the observed 
K' = 17.83 (Ellingson et al. 1996) for lensing implies cB58 has an intrinsic K' ~ 21.59. 32 
of the 74 LBGs for which Shapley et al. (2001) succeed in deriving SED fits have K s within 
0.5 mag of this value. These have a mean stellar mass of 4.0 x 10 10 fyrj M ; a healthy five of 
the 32 have M stars < 3.4 x 10 9 h 2 7 M , but only one has M stars < 0.5 x 10 9 h 2 7 M . Finally, 
following the argument of Pettini et al. (2002), M sta rs ~ 3.4 x 10 9 M would also imply that 
cB58 has already turned a fraction M stars /Mb ary ~ 1/3 of its baryons into stars, roughly 
the level required to account for the enrichment of cB58's ISM to the observed Z = 0AZ & . 
In fact, the stellar/baryonic mass fraction is close enough to this threshold to reinforce the 
argument against previous episodes of star formation; these would tend to have elevated the 
galaxy's oxygen abundance above its observed value. 

Besides charting its past history, we can predict how far into the future cB58 can continue 
forming stars at the present rate before it exhausts its readily available fuel. Combining 
Equations 6 and 8, we have (independent of cosmology) a gas consumption timescale 



For cB58, we find t gas ~ 280^^ Myr, an eminently reasonable value. In nearby spiral 
galaxies, using far-IR luminosity as a proxy for star formation rate, Devereux & Young (1991) 
have showed that the median t gas is ~ 1 Gyr. In contrast, the most luminous local ULIRGs 
have more efficient star formation and t gas approaching ~ 100 Myr (Bryant & Scoville 1999). 
cB58, as befits a starburst of intermediate intensity, falls between these two extremes. 

An important caveat in forecasting the future of star formation in cB58 is that the 
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remaining lifetime of the current episode, tsF(future), ma Y be less than t gas if a large fraction 
of the galaxy's molecular gas reservoir is expelled by an energetic superwind before it can be 
turned into stars. Pettini et al. (2000, 2002) have shown that cB58 is driving a wind in which 
two thirds of the 7 x 10 20 cm -2 neutral column is blueshifted by ~ 170 km s -1 with respect 
to zco (Table 2). The black cores for the strongest absorption lines at this z^ indicate a 
high covering fraction and suggest that much of the outflowing material lies well in front of 
the galaxy's stars along the line of sight. Following Heckman et al. (2000), the mass loss 
rate for a constant-velocity, mass-conserving superwind flowing from some minimum radius 
-R m i n into a solid angle Af2 win d will be 

MqYt- 1 ~ V kpT7 U70 kms-V V7 x 10 20 cm- 2 V \ An J [ ' 

For cB58, setting R min ~ i? uv = 2.0/i^kpc and Af2 wind /47r ~ 0.4 ± 0.2 by analogy with 
local systems (Heckman et al. 2000) yields M W md — (10 ±5) M Q yr _1 , somewhat smaller 
than the current star formation rate. In this simple picture, the superwind should reduce 
^SF(future) below t gas by about a third- an effect whose relative mildness accords with several 
other lines of evidence. Since z^ ^ zco, the bulk of the molecular gas in cB58 cannot 
be participating in the outflow, and the observed enrichment of cB58 implies it has been 
able to retain at least some of the metals produced at earlier times (Pettini et al. 2002). 
Moreover, while tsF(futurc) ~ ^ g as would already fall within a factor of 2-3 of cB58's most 
plausible tsF(past) ~ 100 — 140 Myr, a much lower tsF(futurc) <C tsF(past) would require that we 
have caught the starburst in this system at a special moment near its very end. 

Overall, then, we appear to be observing a galaxy in the midst of converting Mb ary ~ 
10 10 M Q worth of baryons into stars over a period tsF(past) + ^SF(futurc) ^ 420 Myr (see also 
§5.2 below). Both the mass of stars formed (and their age) by the end of this process will 
fall within the ranges inferred from population synthesis analyses of large LBG samples 
(Papovich et al. 2001; Shapley et al. 2001). As expected from considerations of causality 
(Lehnert & Heckman 1996), the duration of the star formation episode will exceed the 
crossing time in cB58, which is roughly t cross ~ 2r co -Da/ cco — 53 \ Myr. By the present 
day, cB58 will be of order a 0.07M* tars system in the local mass function of Cole et al. 
(2001b), 8 with the stars formed in the current burst ranging from (10.9 - 11.3)/i^Gyr 
in age, i.e., indistinguishable from coeval. This projected stellar mass is necessarily quite 
uncertain, given that the superwind may permanently expel some fraction of the system's 
current baryonic mass, while galaxy merging, and ongoing star formation from the atomic 



8 For self-consistency, we assume here that M* tars = 1.4 x 10 11 h 7 M©, as appropriate for the Cole et al. 
(2001b) population synthesis models using a 0.1 — 100 M© Salpeter IMF. 
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gas presumably lying at larger radii, are both likely to occur as time passes. Nevertheless, 
it is suggestive that even this crude estimate lies well within the range of bulge masses 
estimated from .fT-band isophote fitting of local spirals. The samples of Andredakis, Peletier, 
& Balcells (1995) and Moriondo, Giovanardi, & Hunt (1998), for example, scatter (widely) 
about a median Mb u i ge = 3 — 4 x 10 10 M Q , assuming (M/L)b u i g e — 0.6 M Q /L & in K band 
(Moriondo et al. 1998). 

Because cB58 is such a strongly lensed system, it is fair to ask how our conclusions 
depend on the details of the lensing model. The complexity of the S98 model, which is tightly 
constrained by a number of multiply imaged sources in the cluster field, makes a detailed 
examination prohibitive. However, we can at least argue that differential lensing of emission 
at different wavelengths is not likely to have undermined our results. We have already 
assumed without penalty that Mco — M-vv (in §3-1) and that M\jy ~ Mu_ a (above). The 
strong z = correlation between dense gas and dust emission suggests Mfir — Mco is 
also reasonable. More concretely, Shapley et al. (2003) show that cB58 obeys all of the 
same scaling relations (among Ly a equivalent width, interstellar absorption line strengths, 
and interstellar absorption line velocity offsets from Ly a) obeyed by unlensed LBGs. This 
consistency suggests that the fraction (/ arc ~ 0.66) of the background source we are seeing is 
not some unrepresentative slice of the galaxy. In short, differential lensing does not appear 
to have a major impact on the inferred global parameters of cB58, or at least, is not the 
dominant source of uncertainty in our conclusions. It remains possible, of course, to make 
assumptions that will lead to a contradiction of some sort. Adopting a larger conversion 
factor X co c± 7.3 x 10 20 cm~ 2 (Kkms" 1 ) -1 and a smaller age for the current star formation 
episode tsF(past) = 20Myr, for example, the stellar/baryonic mass ratio becomes too small 
(~ 4%) to account for the observed enrichment. However, for the equally (if not more) 
reasonable assumptions made in this paper, no such problems occur; this not only gives us 
confidence that the parameters we have derived are appropriate for cB58, but also hope that 
they may be representative for the LBG population as a whole. 



5.2. Implications for galaxy evolution models 

Our detection of CO (3-2) emission from cB58 provides a set of parameters whose values 
can be compared with those adopted and inferred by models for galaxy evolution at high 
redshift. Basing such comparisons on rather uncertain measurements for a single object 
obviously limits the generality of our conclusions; nevertheless, the exercise is interesting 
as a guide to the analyses that will eventually be possible for larger samples with more 
sensitive millimeter interferometers like ALMA. Our discussion here focuses on comparisons 
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to representative semi-analytic models (Kauffmann & Haehnelt 2000 = KH00; Somerville, 
Faber, & Primack 2001 = SPF01; Cole et al. 2001a = CLBF01) rather than fully numerical 
models, since it is fairly straightforward to compare the assumptions made in the former with 
the relationships among the global properties we have measured for cB58. The first question 
we examine is how well cB58 matches the models' parametrizations for star formation. SPF01 
assume (as in the fiducial ACDM.3 model of Somerville & Primack 1999) that for "constant 
efficiency" quiescent-mode star formation, t gas = 50 Myr at all times. This prescription 
overpredicts the current star formation rate in cB58 by about a factor of five. Other authors 
have chosen to model quiescent-mode star formation using empirical relations seen to hold in 
local disk galaxy samples. We have already noted in §5.1 above that cB58's surface densities 
of (molecular) gas and star formation follow a global Schmidt law; this validates the simple 
assumption made for LBGs by Mo, Mao, & White (1999). Star formation rates also correlate 
with the ratio of gas surface density to a disk's dynamical timescale (Kennicutt 1998; Wong 
& Blitz 2002). For the convention tdyn = 7ri?disk/ , y(-Rdisk), i-e., half of the orbital timescale 
at the disk's outer edge, this relation takes the form locally 

t gas = 31 x 1Q4 /_« x-i /-Rdisk\ MR^yi (12) 



Myr VO.l/ V kpc / V kms 

in terms of a dimensionless efficiency parameter a ~ 0.1 (Kennicutt 1998). For cB58, 
adopting i?disk — 2.0/i^ykpc and v (-Rdisk) — 74 esc? kms -1 gives tdyn — 83 sin i Myr and 
a ~ 0.3 smi \ as observables for the models to reproduce. 

In KH00, the disk dynamical timescale is estimated from the parameters of halos at a 
given redshift under the assumption that disks typically contract by a factor of ten before 
becoming rotationally supported. In this case, 

^L = 9.8xl0»( - X 1 (13) 

Myr Vkms^Mpc- 1 / v ' 

in terms of the Hubble constant as a function of redshift; for our cosmology, H (2.729) = 
282 /10.7 km s -1 Mpc" 1 and tdyn — 35 /i^} Myr. cB58 would then be expected to have a = 
tdynAgas — 0.1 /ioy, compatible with the observations given its uncertain inclination. How- 
ever, a ~ 0.1 also represents almost a factor of four discrepancy with the parametrization of 
quiescent-mode star formation preferred by KH00, which postulates an efficiency declining 
as a(z) oc (1 + z)' 1 . A similar underefficiency is seen in the models of CLBF01, whose 
fiducial case effectively specifies that 

a = 0.005 tt( \ disk - V 5 (14) 
V200kms- 1 / V ; 

for the convention adopted in Equation 12 above and t>di S k ^ ne circular velocity at the half- 
mass radius. Assuming, crudely, that Vdisk — 74 esc i kms" 1 for cB58, Equation 14 predicts 
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an efficiency a ~ 0.0035 esc 15 ? that will only match the observations if % ~ 10°- a scenario 
in which the galaxy's dynamical mass becomes implausibly large. 

Although the above discussion suggests that current semi-analytic models do not ade- 
quately describe cB58 in terms of quiescent-mode star formation, it remains possible that 
a merger-induced burst mode may provide a better account. KH00 model bursts using a 
fixed gas consumption timescale t gas = 100 Myr; this is about a factor of three shorter than 
what we infer for cB58. SPF01 implement a more complex prescription for bursts, in which 
a fraction /consume of the progenitors' original gas mass is turned into stars according to 

SFR (t) = SFR max exp ( - 0.5 (t/t burst ) 2 ) (15) 

over a duration — 2£ DU rst < t < 2tb urs t- Reflecting the results of numerical simulations, 
/consume is set to 0.75 for equal-mass mergers, and to 0.5 for minor (10:1 mass ratio) mergers in 
which neither progenitor has a significant bulge. For cB58, if we fix /consume and assume that 
the current stellar mass within 2.0/i^kpc can be estimated independently of the galaxy's 
star formation history (e.g., as M stars ~ M dyn — M gas ), we can solve for SFR max and tburst as 
well as the age of the current burst. As a major merger in this parametrization, cB58 would be 
observed ~ 14 Myr before the peak of a tburst — 130 Myr burst; as a minor merger, it would be 
observed ~ 35 Myr after the peak of a tburst — 79 Myr burst. SPF01 set tburst oc td yn in their 
models; the td yn — 83 sini \ Myr estimated for cB58 above thus slightly favors a minor 
merger scenario, in which the bulge that could otherwise have suppressed star formation 
(Mihos & Hernquist 1994) might only now be forming (see §5.1). However, since SPF01 
also assume that the quiescent mode of star formation should continue throughout a burst, 
the total star formation rate their model predicts for cB58 is still likely to exceed what is 
observed. 

In the context of the semi-analytic models considered here, then, cB58's star formation 
parameters can be accounted for by a simple Schmidt law; by a quiescent mode whose rate 
depends on the dynamical timescale (provided the efficiency a remains at least as high as seen 
locally; cf. KH00); or by a burst mode triggered by a minor merger (provided rapid quiescent- 
mode star formation can somehow be suppressed; cf. SPF01). Given the observational 
uncertainties, the variety of complex star formation histories that are compatible with cB58's 
rest-UV spectrum (e.g., Leitherer et al. 2001), and the lack of comparison objects, it would 
be premature for us to argue strongly here for or against the use of any of these star formation 
recipes for the LBG population as a whole. cB58 does fulfill the more general expectation- 
noted by both KH00 and SPF01- that star-forming galaxies at high redshift should be 
relatively gas-rich. SPF01 in particular predict that LBGs should have high mass fractions 
of cold (molecular) gas, with the median / gas ~ 0.5 for LBGs with < 25.5. For cB58, 
whose observed V magnitude (Ellingson et al. 1996) corrected for lensing implies V ~ 24.4, 
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the most likely / gas ~ 0.66 falls well inside the envelope of the authors' mock Hubble Deep 
Field catalogue (as is also the case for several of the galaxy's other parameters; see their 
Figure 16). The KH00 models predict similarly that / gas for major mergers increases to 
~ 0.6 by z ~ 2.7, although the evolution will not be as strong for a ACDM cosmology. 

A second question our observations of cB58 help address relates to the uncertain mass 
scale of LBGs. Direct estimates of stellar masses from optical/near-IR photometry (Papovich 
et al. 2001; Shapley et al. 2001) and dynamical masses from nebular emission line widths 
(Pettini et al. 1998, 2001) reveal Mb ary ~ 1O 1O M concentrations in the centers of LBGs, 
as anticipated by early analyses of LBG sizes and luminosities alone (e.g., Lowenthal et al. 
1997). In this context, it is not a surprise that we estimate Mb ary ~ 10 10 M in cB58 from 
observations that (by assumption- most explicitly in our invocation of Aico — -Muv) again 
trace only a galaxy's bright core. The fact that the CO (3-2) linewidth is within 10% of 
the Ha linewidth measured by Teplitz et al. (2000)- <xco = 74kms~ 1 vs. ana = 81kms _1 - 
does supply useful evidence that Ha is not unduly affected by the superwind, although it 
remains possible that neither Ha nor CO emission traces cB58's full rotation curve. The 
more unique contribution here comes from our estimate of a total star formation timescale 
(^SF(past) + ^SF(futurc) ^ 420 Myr for a constant star formation rate; 4 x tburst ~ 320 Myr for a 
merger-triggered burst a la SPF01) that is a substantial fraction of the time elapsed over the 
2.5 < z < 3.5 range probed by the Lyman break selection function (t^z ~ 800 Myr for our 
cosmology). Given that the star formation timescale can be extended if the gas reservoir is 
replenished by inflow, recycling (Kennicutt, Tamblyn, & Congdon 1994), or future merging, 
we can conclude- provided cB58 is representative and does not suffer substantial mass loss- 
that the duty cycle for star formation in LBGs at this redshift will be of order 50%. 

The strong clustering of LBGs implies that they are a highly biased galaxy population 
(Giavalisco et al. 1998; Adelberger et al. 1998; Arnouts et al. 1999; Giavalisco & Dickinson 
2001; Porciani & Giavalisco 2002; Foucaud et al. 2003), and thus that they must occupy 
the most massive halos at z ~ 3. Exactly how to associate LBGs with halos of mass M- 
i.e., the precise form of the halo occupation function N g (M)- remains an open issue. If 
N g = 1 is assumed, a fiducial M can be estimated for a given LBG sample by identifying the 
mass scale at which halos have the same number density and bias. Adelberger et al. (1998) 
and Giavalisco & Dickinson (2001) use this approach to infer M ~ 1.1 x 10 12 for the 
halos hosting LBGs with 1Z < 25.5 (like cB58 after correction for lensing). The associated 
baryonic masses should then be M bary ~ 1.5 x 10 11 M Q , 9 which even if reduced to match 
lower estimated clustering strengths (Adelberger et al. 2003) will be in excess of the direct 



9 We assume M/M hmy = Q M /^b for Q h = OMh^ (Buries & Tytler 1998). 
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measurements for most individual LBGs by at least a factor of a few. Packing multiple 
LBGs into a single massive halo would in principle reduce this discrepancy. However, while 
the high close pair fraction at z ~ 3 does favor dN g (M)/dM > (Wechsler et al. 2001; 
Bullock, Wechsler, & Somerville 2002), the apparent dependence of clustering strength on 
number density (Giavalisco & Dickinson 2001; Foucaud et al. 2003) argues strongly that 
in fact dN g (M)/dM < (Bullock et al. 2002). Our estimate of a ~ 50% duty cycle for 
star formation further suggests that the remaining baryons in a halo hosting one LBG also 
do not lie in some large population of cohabiting, "dormant" LBGs in which starbursts 
have faded or not yet begun. This consideration further disfavors scenarios in which star 
formation episodes in an individual LBG are extremely brief (e.g., triggered by frequent 
sub-halo collisions: Kolatt et al. 1999). Instead, the remaining baryons in a given LBG's 
host halo are likely to be found in the galaxy's own lower surface brightness outskirts, in 
cohabiting galaxies of significantly lower mass, or outside galaxies altogether (where they 
may have been banished by a vigorous superwind). 



5.3. Millimeter line/continuum ratios at z ~ 3 

After SMMJ1401 1+0252 (Frayer et al. 1999; Ivison et al. 2000), cB58 is only the second 
galaxy without an active nucleus to have been detected in CO emission at high redshift. 
Among all detected sources, cB58 is also fairly unusual in not having been prioritized for 
CO observations on the basis of a previous measurement of strong dust emission in the 
(sub)millimeter; the radio galaxy 53W002 (Scoville et al. 1997) is perhaps the only other 
member of this category. Because it does not reflect common selection biases, then, it is 
of interest to compare cB58's line/continuum ratio to those of other systems at comparable 
redshifts. Due to the large uncertainties in dust temperature and X co conversion factor, this 
ratio may only poorly translate to a star formation efficiency or a gas/dust ratio. However, it 
is still interesting from a practical standpoint to see how much source-to-source variation in 
global ISM parameters exists at high redshift, since this will drive strategies for deep surveys 
with the next generation of millimeter interferometers. 

In order to limit the impact of CO excitation uncertainties on our comparison, we will 
restrict our comparison to sources in the range 2.3 < z < 2.8 that have been observed in 
the CO (3-2) line and in continuum at 1.2 mm. The observed 1.2 mm flux density (corre- 
sponding to rest wavelengths 316 — 363 /im over this redshift range) is a generally reliable 
(albeit somewhat redshift-dependent) measure of thermal dust emission. 10 cB58 has a ra- 



'We exclude the exception to this rule, MG 0414-0534, whose extremely strong 1.2 mm flux density is 
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tio = Fqo/Si.2 ^ 350 km s -1 ; if we assume as discussed in §4.2 that 0.22 mJy of its 

apparent 1.2 mm flux density actually comes from the foreground cD, then 7^l/c becomes 
~ 440kms _1 . The seven comparison objects in the compilation of Beelen et al. (2003) 11 
have values of 72-l/c ranging from 300 km s -1 (J 1408+5628: Omont et al. 2003; Beelen et 
al. 2003) at the low end to 700 km s" 1 (53W002: Alloin, Barvainis, & Guilloteau 2000) and 
800 kms" 1 (SMM J14011+0252) at the high end. With cB58 occupying the opposite end of 
this interval from the only other galaxy not hosting an energetically dominant active nucleus 
(SMM J14011+0252), and the only other galaxy not selected on the basis of its dust emis- 
sion (53W002), there is no evidence that either selection criterion significantly biases the 
line/continuum ratio. It is in fact rather remarkable that all eight systems have 72-l/c values 
agreeing to within better than a factor of three. Two quasars with CO (3-2) upper limits in 
the survey of Guilloteau et al. (1999) would in principle correspond to 7^l/c < 300 km s -1 ; 
however, since quasar nondetections often result from uncertain CO redshifts rather than 
weak CO lines, these may not be the outliers they appear. 

6. Conclusions 

From our observations of the gravitationally lensed Lyman break galaxy cB58, we esti- 
mate that within a radius of 2.0 /i^kpc, it most likely has 

1. an age for the current star formation episode tsF(past) ~ 140 Myr; 

2. a current star formation rate = (23.9 ± 0.7) M Q yr -1 ; 

3. a formed stellar mass M sta rs ~ 3.4 x 10 9 M ; 

4. a molecular gas mass M gas = 6. 6^4:3 x 10 9 M ; 

5. a lifetime for continued star formation tsF(futurc) < ^ g as = 280^ 2 g°,Myr, lower to the 
extent that molecular gas is expelled by the observed superwind before it can be turned 
into stars; 

6. a molecular gas mass surface density S gas = 540^35°, M Q pc~ 2 ; 

7. a dynamical mass M dyn = l-Ol 9 ;^ x 10 10 csc 2 i h Q \ M Q ; and 



dominated by nonthermal emission (Barvainis et al. 1998). 

11 F10214+4724, the Cloverleaf, 53W002, SMM J02399-0136, Q1230+1627B, SMM J1401 1+0252, and 
J 1408+5628. 
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8. a gas mass fraction / gas = 0. 661^52 sm ^ ^o.7- 

We judge that these estimates should not have been severely compromised by differential 
lensing, or by the possible effects of nonthermal continuum or (lensed) CO line emission from 
a nearby source X. We conclude that cB58 is partway through an episode of star formation 
lasting several dynamical timescales, and in the absence of substantial mass loss, gas inflow, 
or further merging will evolve into az = system whose mass and coevality are those of a 
spiral bulge. 

To the extent that cB58 is a typical LBG apart from its magnification, its global pa- 
rameters place new constraints on the permissible recipes for star formation in semi-analytic 
models of galaxy evolution. Its current star formation rate, molecular gas mass, and inferred 
stellar mass are at least reasonably compatible with the expectations for a Schmidt law, for a 
quiescent mode of star formation tied to the local dynamical timescale, and for a burst mode 
of star formation (here triggered by a minor merger), although no published model offers 
perfect agreement. Its high gas mass fraction agrees with theoretical predictions that / gas 
should increase with redshift. The duration of cB58's current star formation episode appears 
to be long relative to the redshift range probed by the Lyman break technique at z ~ 3; this 
argues for a high duty cycle for star formation, and against scenarios in which the halo host- 
ing an observed LBG is over time "lit up" by many different objects undergoing brief bursts. 
Finally, from a comparison with other z ~ 2.5 systems detected in CO (3-2) and dust emis- 
sion, we arrive at the useful empirical conclusion that the millimeter line/continuum ratio is 
not a strong function of nuclear activity or dust luminosity. 

We are grateful to IRAM staff members at the PdBI and in Grenoble, particularly 
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as for her insights on lens models, and acknowledge helpful comments from an anonymous 
referee and useful discussions with Andrea Cimatti, Pierre Cox, Bruce Elmegreen, David 
Hogg, Guinevere Kauffmann, Masami Ouchi, Max Pettini, Alice Shapley, Rachel Somerville, 
and Harry Teplitz. This work has made use of the NASA/IPAC Extragalactic Database 
(NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with NASA. 
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A. Revisiting the S98 lensing model 

S98 model the total deflection potential of the lens as a sum of elliptical potentials due to 
individual galaxies and the cluster as a whole. The ith term in this sum is centered at some 
position (xi,yi) and is characterized by a normalization Q Ei (recoverable from a tabulated 
velocity dispersion <7j according to their equation A3), a scale length Q, and an ellipticity ef 

<j> t (x, y) = 6 E , (C 2 + ^ (x - x t f + I±* (y - ytf) * (Al) 

In gravlens, a softened power-law potential can be specified in terms of an index a iy normal- 
ization bi, scale length Sj, and ellipticity ef. 

& (x, y) = h (s 2 + (x- Xi f + ^ (y - y,) 2 ) 2< *° (A2) 

Equality between Equations Al and A2 holds for ctj = 1 and 



6* = ©^(^K (A3) 
1 + e A § 



(A5) 



2e, ; 



1 + ei 

Our version of the S98 model is constructed in an (x, y) coordinate frame with axes parallel 
to the rows and columns of the WFPC2 chip and an origin at the position of the cluster 
cD. We have measured lens-plane galaxy positions {x i: yi} from the F675W image of S98, 
used Equations A3 through A5 to rescale the {6e,i, d, values specified for their model 
MVIa, and rotated position angles for the cluster and cD major axes to match our coordinate 
convention. To confirm that the resulting gravlens version is a reasonable facsimile of the 
original model, we have examined the lensing of a background point source which would 
appear at the observed position of image A2. Our model specifies that the image of a point 
source at A2 has undergone a factor of 2.38 magnification, in excellent agreement with the 
factor of 2.39 estimated by S98. Our model also predicts that two additional images of the 
same background source should appear at the locations marked by filled circles in Figure 
5, i.e., within the cB58 arc. Exact reproduction of the arc produced by an intrinsically 
extended (rather than point) background source is beyond the scope of this paper. 

We have also checked that the projected mass density specified by our reparametrization 
of the S98 lensing model can roughly account for the positions and magnifications of the 
additional (B, C, and W) sets of multiple images in this field. Because the corresponding 
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background sources do not necessarily lie at the same z ~ 2.73 as cB58, it is necessary 
to apply a further rescaling to the deflection potential for each one. Conveniently, the 
cosmological dependence of the deflection potential can be factored out of the normalizations 
for the individual terms. For source and lens redshifts z a and z\ and the critical surface density 
for lensing E crit , 

, ( x v -i 4ttG D a (z h 0) D A (z s ,0) 

0, (x, y) oc S cr = — (A6) 

where we write D A (z2, z±) to denote the angular diameter distance to a source at redshift 
z 2 as seen from a lower redshift z 1 . The correct rescaling for the jth background source at 
redshift Zj is thus to map b; L b; L dj for each i, where dj is a lensing strength defined as 

, f v Da(z v 0.37) r D A (2.73, 0.37) ] -i 

d 3 ( Zj , n M , n A ) = Da( ^ q) [ Da(2 , 73;0 ) J ( A? ) 

The precise mapping between and Zj depends on the choice of cosmology. Since the S98 
model MVIa solves for the {dj} explicitly (corresponding, in our cosmology, to zb — 3.05, 
zc — 3.57, and zw — 2.85), we can rescale our gravlens deflection potential directly. The 
derived image positions shown in Figure 5 are in generally good agreement with the data, 
although gravlens predicts that there should be one additional image each of both the C and 
the W background sources (only the latter is predicted by S98). Note that the assumption 
of z x = 1-48 (see below) fixes dx = 0.874 in our cosmology. 

In exploring the possibility that source X could be gravitationally lensed CO (2-1) emis- 
sion from a source at z x ^ 1.48 or CO(3-2) emission from a source (other than cB58) at 
z x — 2.73 (see §4.1 above), the strongest constraint is imposed by the absence of counterim- 
ages within the 50% response of the PdBI primary beam at 92.96 GHz (see, e.g., the lower 
panel of Figure 6). Indeed, the second-strongest feature in the entire map is a 4<r peak lying 
~ 57" to the southeast of X (not shown), where the primary beam response is only ~ 3%. If 
this were a true counterimage, its implied flux density of 20.4 ±5.0 (±3.0) mJy would require 
that it be more highly magnified than X by a factor between 15 and 36- an unlikely situation 
at a position where the lens model predicts only a single image and a magnification ~ 1.1 
for either choice of z X - The counterimages of X that the model does predict lie at much 
smaller radii. For both zx — 1-48 and zx — 2.73, we have used gravlens to predict the image 
positions and magnifications that result if we place X at each of the nine grid points defined 
by (a ± Act, 5 ± AS) in terms of the astrometric uncertainties listed in Table 1. Typically, 
we find X has one or more counterimages that are at least as strongly magnified, far enough 
away not to be confused with X, and close enough to our phase center that they would hardly 
suffer any primary beam attenuation. Such counterimages would be detected if they existed 
at all; their absence instantly invalidates the particular models in which they appear. 
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In a few cases, the lensing geometry can conspire to place X and a bright counterimage 
X2 close to each other on opposite sides of a critical curve. The upper panels of Figure 6 
show the most promising examples of such configurations. For Zx — 1-48 (at left), placing X 
at (a — Aa, 5) yields a net magnification of X+X2 three times higher than that of the next 
brightest counterimage in the field. For zx — 2.73 (at right), the magnification of X+X2 
that results from placing X at (a + Aa, 5) is a factor of twelve higher than that of the next 
brightest counterimage- entirely consistent with the manifestation of only a "single" source 
in the lower panel of Figure 6. Exactly how realistic such models are depends strongly on the 
precise shape of the lens mass profile at small radius. With Kassiola & Kovner (1993) having 
shown that elliptical potentials translate to unphysical, peanut-shaped isodensity contours 
for ellipticities e > 0.18 (using the gravlens definition), for example, our use of e = 0.28 for 
the cluster and e = 0.23 for the cD suggests there are reasonable grounds for caution. 
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Table I. Positions and positional uncertainties. 





^obs 


a 


Aa 


5 


A5 




Source 


(observed) 


(J2000) 


(J2000) 


(J2000) 


(J2000) 


Reference (s) 


phase center 




15:14:22.25 




±36:36:24.40 






cB58 a 


CO(3-2) 


15:14:22.218 


±0.050 


±36:36:24.79 


±0.39 


1 




optical 


15:14:22.271 e 


±0.010 


+36:36:25.2 e 


±0.12 


2,3 


cD fe 


6.1 cm 


15:14:22.47 


±0.002 


±36:36:20.8 


±0.02 


2 




21 cm 


15:14:22.443 


±0.075 


±36:36:20.68 


±0.90 


4 


S c 


optical 


15:14:22.578 e 


±0.010 


±36:36:16.8 e 


±0.12 


2,3 


X 


3.2 mm 


15:14:22.455 


±0.049 


±36:36:20.90 


±0.38 


1 



a cB58 = MS 1512.4±3647: PPP 101120 in Abraham et al. (1998). 

b cD = MS 1512.4±3647: PPP 101094 in Abraham et al. (1998). 

C S (S98) = MS 1512.4±3647: PPP 101068 in Abraham et al. (1998). 

c Calculated by identifying the Laurent-Muehleisen et al. (1997) 6 cm radio core with the cD, 
and using the Abraham et al. (1998) offsets from the optical cD position. 

References. — (1) this paper (2) Laurent-Muehleisen et al. (1997) (3) Abraham et al. (1998) 
(4) Becker et al. (1995) 
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Table 2. Measured redshifts for cB58, from red to blue. 



Quantity 


Rest feature (s) 




Redshift 


Velocity 
(kms -1 ) 


Reference 


z b 

stars 
stars 

zco 

~c 

^abs 


mean of 10 (AA 3726 - 

CIII A 2297 
mean of 10 (AA 1294 - 

CO (3-2) 
mean of 35 (A A 1134- 


6583) 
1718) 
2577) 


2.7290 ± 0.0007 

9 797«+0.0001 
z -' z ' u -0.0003 

2.7268 ± 0.0008 

9 79«c;+0-0004 
z,. i z,uo_ 0005 

2.7244 




-110 
-180 
-200 
-370 


1 

2 
3 
4 
2 



a Velocity relative to zun in the source rest frame. 

b We list both reported estimates of z stars , although Pettini et al. (2002) express a 
preference for -2 stars = 2.7276. 

c Redshift of deepest interstellar absorption. Pettini et al. (2002) note that absorption 
at some level occurs over a velocity range of 1000 km s -1 . 



References. - (1) Teplitz et al. (2000) (2) Pettini et al. (2002) (3) Pettini et al. 
(2000) (4) this paper 
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10 -10 
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Fig. 1. — Zeroth moment of the core of the CO(3-2) line in cB58. Coordinates are relative to 
the PdBI pointing center listed in Table 1. Emission has been averaged (without blanking) 
over a 60 MHz (~ 194.0 km s -1 ) window centered at 92.789 GHz (185.9 km s _1 blueward of 
zhii), i.e., over four of the 48.5 km s -1 channels in Figure 2 below. Contours are multiples 
of 0.3 mJy beam -1 for a synthesized beam of 8.2" x 4.8" at P. A. 73° (shown at lower left); 
negative contours are dashed. Interpreted as line emission, the 1.7± 0.3 mJy beam -1 peak at 
the phase center corresponds to a CO(3-2) flux of 0.33 ± 0.06 Jy kms" 1 . No primary beam 
correction has been applied. 
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Fig. 2. — CO(3-2) spectrum of cB58 with a Gaussian fit overlaid. Velocity channels are 
independent and 48.5 km s -1 wide; the zero point of the velocity axis corresponds to zhii = 
2.7290 derived by Teplitz et al. (2000). Gaussian a co = FWHM co /2.355 = 74 ± 18km s -1 . 
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Fig. 3. — Pseudo-continuum images at three different 3 mm frequencies. Synthesized beams 
are plotted at lower left; crosses show ±la positional uncertainties for the centroid of CO(3- 
2) emission from cB58. Contours are multiples of 0.24 mJy beam -1 , which corresponds to 
la for the left and right panels and 1.5a for the center panel; grey contours are negative. 
The 5.4a source seen to the southeast of cB58 in the center panel is called X in the text. No 
primary beam correction has been applied. 
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Fig. 4. — 3 mm "spectrum" of X with the same velocity channels as Figure 2. Dashed lines 
show the approximate windows averaged to produce the center and left panels of Figure 3. 
Solid line at bottom shows the velocity range spanned by the core of cB58's CO(3-2) line as 
used to plot Figure 1. 
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Fig. 5. — The S98 lensing model as implemented in gravlens, with coordinate offsets relative 
to the location of the cD and parallel to the rows and columns of the WFPC2 chip. For each 
image set (A = circles, B = squares, C = diamonds, W = triangles), open symbols show 
observed images and solid symbols show predicted images. One model image in each set 
(A2, Bl, CI, and Wl) has been forced to agree with its observed location; the reasonable 
agreement between model and data for the other images of the corresponding background 
source verifies that the model works reasonably well. The thick lines show the image-plane 
critical curves for a background source at z ~ 2.73 (compare e.g., Figure 14 of S98), only 
directly relevant for the A image set that includes cB58; the thin circle shows the 50% 
response of the PdBI primary beam at 92.96 GHz, whose diameter is ~ 52". 
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Fig. 6. — Lower panel: same as center panel of Figure 3, but with a larger field of view and 
contours that are multiples of la (0.16 mJy beam^ 1 ). Upper panels: possible lens models 
for z x — 1-48 (left) and zx — 2.73 (right), in coordinates defined relative to the cD position 
and rotated into alignment with the WFPC2 chip. Small circles at origin show the nominal 
position of X; filled circles show the positions of predicted images and are labelled with 
magnification factors. The thick lines show the image-plane critical curves for a source at 
zx- The thin circle in all three panels shows the 50% response of the PdBI primary beam 
at 92.96 GHz. 



